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E. coli bacteria swim following a run and tumble pattern. In the run state all flagella join in a 
single helical bundle that propels the cell body along approximately straight paths. When one or 
more flagellar motors reverse direction the bundle unwinds and the cell randomizes its orientation. 
This basic picture represents an idealization of a much more complex dynamical problem. Although 
it has been shown that bundle formation can occur at either pole of the cell, it is still unclear whether 
this two run states correspond to asymmetric propulsion features. Using holographic microscopy 
we record the 3D motions of individual bacteria swimming in optical traps. We find that most cells 
possess two run states characterised by different propulsion forces, total torque and bundle confor¬ 
mations. We analyse the statistical properties of bundle reversal and compare the hydrodynamic 
features of forward and backward running states. Our method is naturally multi-particle and opens 
up the way towards controlled hydrodynamic studies of interacting swimming cells. 


INTRODUCTION 

Escherichia coli is a peritrichously flagellated rod¬ 
shaped bacterium with no apparent signs of polarity. It 
swims at about 20 jam/s following a characteristic run 
and tumble pattern [T]. When all flagella rotate in the 
same direction they join together near a cell pole forming 
a compact helical bundle that propels the cell body along 
approximately straight runs. The bundle, however, has 
a finite lifetime of about 1 s after which it disentangles 
due to one or more flagellar motor changing direction [2] . 
In this transient state, called tumble, the uncoordinated 
motions of flagella cause a rapid and random reorienta¬ 
tion of the cell body axis. After about 0.1 s all flagella 
resume to the same rotation direction and a new bundle 
forms. In a perfectly symmetric cell, bundle formation 
should occur with equal probability at either cell pole 
and result in equal propulsion forces. Conversely, polar 
flagellated bacteria like P. putida have been observed to 
swim in opposite directions with very different speeds [3] . 
However, even in the case of an apparently symmetric cell 
like E. coli , division by binary fission allows to distin¬ 
guish between a new and an old pole. An age asymmetry 
between the two poles will naturally give rise to a corre¬ 
sponding asymmetry in all pole properties that depend 
on the cell age. For example, the chemoreceptor Tsr has 
been shown to be more abundant at the old pole [4] . An 
accurate analysis of flagella distribution over the cell sur¬ 
face also shows an increased density of flagella over the 
old pole cell half [5]. In addition, since the average num¬ 
ber of flagella is typically only about 6 [5], large devia¬ 
tions from a symmetric flagellar distribution are expected 
to occur just by random chance. These observations lead 
naturally to the question whether the motility properties 
of E. coli display any corresponding asymmetries, what 
is their nature and entity and if there’s any particular ad¬ 


vantage that an asymmetric swimming could offer. One 
of the first polar features that was recognised in the run 
and tumble motion of E. coli was that a large number of 
cells show a marked tendency of forming the bundle at 
one specific pole after a tumble [6]. Later work showed 
that the old pole is usually the preferred pole for bundle 
formation and attributed this to a larger average num¬ 
ber of flagella over the old pole cell half [5]. Long term 
observation of tumbling E. coli in a double optical trap 
confirmed the presence of a preferred pole for bundle for¬ 
mation in the majority of cells [ 7 ]. The latter work also 
reported the presence of two distinct body rotation fre¬ 
quencies corresponding to a bundle forming at one pole 
or the other. However, the presence of two corresponding 
distinct propelling forces has never been reported so far. 

A direct measurement of the propelling force is feasi¬ 
ble in principle using calibrated optical traps [8]. Using a 
two trap setup, combined with back focal plane detection 
of cross-polarised transmitted beams, is a common and 
highly accurate way for measuring forces between two 
trapped spherical objects [9]. Although a two trap con¬ 
figuration has been successfully applied for trapping and 
alignment of rod-shaped bacteria [7] , force measurements 
are not straightforward when trapping a single elongated 
body. An alternative approach consists in trapping a cell 
near the coverslip with a single beam while simultane¬ 
ously applying an external flow to align the cell on the 
focal plane and gauge propelling forces with viscous drag 
m- The presence of the flow however could strongly bias 
bundle formation towards one pole. Moreover, due to the 
proximity of a solid wall, hydrodynamic wall effects can 
strongly affect force measurements and swimming effi¬ 
ciency. 

Here we demonstrate that a combination of single 
beam optical trapping and 3D holographic microscopy 
provides an ideal tool to investigate in the bulk fluid 
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the propulsion features of the two possible run states 
together with a careful examination of the statistics of 
transitions between them. We find that the majority of 
cells displays two distinct run states depending on which 
side the bundle forms. The average difference between 
the two corresponding propelling forces is 35% with the 
largest deviations exceeding 100%. 


METHODS 
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FIG. 1. (a) Optical setup: three RGB LED sources are used 
for 3-axes DHM volumetric reconstructions of individually 
trapped cells. The trapping IR laser (shown in gray) is phase 
modulated by a SLM and sent to the entrance pupil of a 
lOOX, NA 1.45 oil immersion objective by a beam splitter 
(BS). A second 532 nm laser is coupled to the objective and 
used to stop bacterial motility by irreversible photodamage, 
(b) A typical RGB fringe pattern produced by a trapped cell 
located 7 /xm above the focal plane, (c) Gorresponding 3D 
reconstruction of the cell body. 


E. coli cells (MG1655) were streaked on a Petri dish 
containing 1.5% of agar and tryptone broth (1% tryp- 
tone and 0.5% NaCl). A single-colony was inoculated 
into tryptone broth and grown overnight at 32°C. The 
saturated culture was then diluted 1:100 (50 jal in 5 ml) 
into fresh medium and grown for 4 hours at 32°C with 200 
rpm rotation. Bacterial cells were then harvested from 
culture media by centrifugation at 2200 rpm for 10 min¬ 
utes at room temperature. The pellet was resuspended 
by gently mixing in trapping buffer composed of 100 mM 
Tris-HCl, 70 mM NaCl, 2% (wt/vol) glucose and an oxy¬ 
gen scavenging system [7]. The cells were washed three 
times to replace growth medium with trapping buffer. 
Trapping buffer sustain bacteria motility and reduce the 
oxidative damage to the cell induced by trapping. Indi¬ 
vidual cells were trapped using an infrared (A=1064 nm) 
holographic trapping system m based on a LCOS Spa¬ 


tial Light Modulator (Hamamatsu X10468) Figure [^a). 
Using an oxygen scavenging system trapped bacteria can 
display a high motility for hours. 3D volumetric recon¬ 
structions of trapped bacteria are obtained using a 3- 
axes/3-colours modification of in-line Digital Holographic 
Microscopy (DHM) [12]. In conventional in-line DHM 
the sample volume is reconstructed by numerically back- 
propagating the complex field encoded in the interference 
pattern between scattered and incident light [T3| . In our 
setup we replaced the single on-axis illumination beam 
with three LED collimated light beams of different wave¬ 
lengths (Prizmatix UHP-Mic-LED-630, UHP-Mic-LED- 
520, UHP-Mic-LED-450) chosen to match the spectral 
response of the three colour channels of a RGB colour 
camera (Easier avAlOOO-lOOgm). The sample scatters 
the light of the three beams producing three indepen¬ 
dent fringe patterns. The field intensities corresponding 
to the three colour channels are reconstructed in 3D and 
then are overlapped to obtain a volumetric image. Fig- 
ure[^b-c) shows a RGB hologram with the corresponding 
3D reconstruction of a trapped E. coli . As in standard 
DHM, good numerical reconstructions are obtained when 
the scattering object is located a few microns above the 
focal plane [H HI- To meet this requirement, we use 
the SLM to move the trapped cells 7 /xm above the focal 
plane. 


RESULTS 

We trapped individually 53 cells and observed their 3D 
motion at 100 frames per second. Total observation time 
for each swimming cell was 5 minutes during which no 
motility reduction was observed. Cells were trapped 40 
/xm above the cover slip so that hydrodynamic coupling 
to the glass surface was negligible. 

Rod shaped objects like E. coli spontaneously align 
along the beam axis maximising the overlap between the 
prolate cell body and the elongated focal distribution of 
light. In addition to trapping gradient forces, that oppose 
cell displacements away from the trap center, gradient 
torques arise restoring the cell axis onto the beam axis. 

During the run state, a flagellar bundle is formed on 
one pole pushing the cell body away from the trap equi¬ 
librium axial position z = 0. Due to a slight misalign¬ 
ment between the bundle and body axes, the cell body 
precesses around 2 :. The time evolution of cell position 
and orientation were extracted from the first and second 
spatial moment of volumetric images. Figure |^b) shows 
typical time traces obtained for the axial components of 
cell body position z, angular velocity and polar an¬ 
gle 0 of the body axis. Synchronous transitions between 
two states are clearly visible in the 2 : and traces cor¬ 
responding to the two possible run modes in an optical 
trap: i) bundle pointing to the negative ^ direction and 
cell swimming upwards with positive 2 : components for 
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FIG. 2. (a) A trapped cell can form a bundle on either of its 
two poles. The rotating bundle at one pole pushes the cell 
body towards the opposite pole. The two swimming states 
also give rise to opposite cell body rotations around z. (b) 
Time traces of the axial cell body displacement from the trap 
equilibrium position z, rotational frequency and cell body 
polar angle 0. (c) 2D histogram obtained from the z, traces 
shown in (b). 


vation time window. We counted a total of about 5400 
tumble events, of which 34% led to a bundle reversal. 
This value is in between the previously reported figures 
of 56% [6] and 16% [ 7 ]. For the particular cell shown in 
Figi we counted 99 tumble events of which 34 lead to 
a bundle reversal. The cell spends about the same time 
in the two run states with an average run duration of 3 
s. Fig. |3] shows the cumulative run length distributions 
for the two states of 4 cells. Run length distributions 
display a high variability among different cells (different 
colours) but also between the two running states of an 
individual cell (solid and dashed lines). The distribution 
of running times is exponential for run lengths shorter 
than 5 s but shows a systematic excess of long runs, as 
also found in mm- In particular, our data seem to show 
that the fat tail at long run lengths is more pronounced 
in one of the two running states. Trap calibration is per- 



displacement and angular velocity, ii) bundle pointing to 
the positive z direction and cell swimming downwards 
with negative z components for displacement and angu¬ 
lar velocity. A closer look to a zoomed 10 s region also 
reveals the presence of tumble events that do not lead 
to a bundle reversal (dashed lines in Fig. |^b)). These 
events can be discriminated from bundle reversals (solid 
red lines in Fig. ib)) since both 2 ; and 11^, after a ^0.3 
s transient where large fluctuations around zero are visi¬ 
ble, both return to their initial value rather than chang¬ 
ing sign. These transient tumbling states however are 
not visible in the 2 ; and histogram (Figure |^c)) for 
two reasons: i) the systems spends a negligible amount 
of time in a tumble compared to to the time spent in the 
two run states, ii) the values of 2 : and are not well 
defined during a tumble. The time trace of the cell polar 
angle 0 reveals the presence of reorienting torques acting 
on the cell body during tumble and reversal events (Fig. 
[^c)). However, the restoring torque applied by the op¬ 
tical trap always prevents the cells from reorienting and 
constraints the polar angle below 5 degrees. 39 bacteria 
out of 53 displayed bundle reversals within our obser- 


FIG. 3. Gumulative run length distributions. Different 
colours represent different cells. Solid and dashed lines dis¬ 
tinguish between the two possible run states corresponding 
to having the ffagellar bundle at each of the two cell poles. 
Black line is the logarithmic average of all the reported curves 
whose corresponding exponential ht is shown in gray. 

formed independently for each cell after a short exposure 
to a focused 532 nm beam that causes photodamage and 
stops motility [EHII]. Using our 3D imaging technique 
we can track the Brownian fluctuations of the cells posi¬ 
tion and orientation and reconstruct the optical potential 
assuming Boltzmann statistics (see Appendix). Once the 
trap is calibrated we can directly convert z displacements 
into propelling forces. In particular we want to address 
the open question: is there any change in the propelling 
force associated with rebundling of flagella on the oppo¬ 
site cell pole? Before we go further with data analysis it 
is useful to consider a simple hydrodynamic scheme for 
the swimming problem. Following [18] we model a swim¬ 
ming cell as the combination of two rigid units: the cell 
body and the flagellar bundle. The two units are rigidly 
connected but they are free to rotate around a common 
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axis. We can formulate the corresponding hydrodynamic 
equations in terms of a resistance problem connecting 
force and torque components along the swimming direc¬ 
tion to the corresponding velocity and angular frequency 
components: 


pb = A^U 

( 1 ) 


( 2 ) 

pS = A^U + 

( 3 ) 

Tf = B^U + D^n^ 

( 4 ) 


where the superscripts b and / refer respectively to cell 
body and flagellar bundle. The symbol A represents the 
translational drag coefficient connecting the linear speed 
U to the force F that is needed to move a body in a 
purely translatory motion. Similarly D connects angular 
speed to torque in pure rotations. For bodies that, like 
the helical bundle, do not posses three mutually orthog¬ 
onal planes of symmetry, a third coupling coefficient B is 
required connecting the angular speed to the force that 
is required to prevent them from translating. 

For a trapped cell the swimming speed U vanishes 
while the total external force is provided by the optical 
trap: 

pb + F^ = Ftrap = Bfnf (5) 

In other words the trapping force balances the thrust 
Fthrust = from the bundle and prevents the cell 

from swimming. We measure an average thrust of 0.14 
pN which is lower than the previously reported values 
0.57 pN and 0.4 pN respectively in [10] and [T9|. This is 
consistent with the lower average swimming speed that 
we measured on our samples 10 /im/s) compared to 
an almost double value 20 /im/s) reported in [TOlfT^ . 
We believe this discrepancy could be due to differences 
in strain, buffer composition and sampling of bacterial 
tracks. 

In order to quantify the degree of variability of a 
generic positive quantity Q between the two run states, 
denoted by the subscript + and —, we use the normalised 
difference A[(5] = 2(|(5+ — (5-|/(Q+ + Q-)) where an¬ 
gle brackets denote averages over different cells. Fig. 
I^b) reports a histogram of the normalised force differ¬ 
ences between the two states. For the thrust we obtain 
A [Fthrust] =0.35 that is much larger than the average 
force fluctuations between different runs in the same di¬ 
rection (10%). In other words, run and tumble dynam¬ 
ics in E. coli is characterised by two different and well 
defined run states corresponding to flagella bundling at 
one of the two cell poles and resulting in propelling force 
differences that are on average 35% with tails extend¬ 
ing over 100%. This is in contrast with what found in 
another peritrichously flagellated bacterium B. subtilis 
where very similar swimming speeds have been measured 


m- The more symmetric swimming behavior of B. sub¬ 
tilis could be related to a larger number of flagella 30). 
As shown in Fig. [^c), large differences are also observed 
between the torque values in the two run states. Since 
does not change upon bundle reversal we have, using 
([^, A[T^] = A[f2^] = 0.45 that is again much larger than 
torque fluctuations between different runs pointing in the 
same direction (15%). It is interesting to note that, using 
the torque free condition = —T^ and applying A[-] to 
both sides of we get A[T^] = A[T^] = A[F^]-h A[f]^]. 
If we assume, as reported in [ 7 ], that the bundle rota¬ 
tion frequency does not change much upon reversal then 
A[(2-^] ~ 0 and therefore we obtain A[T^] ^ A[D^. This 
seems to suggest that the observed torque differences be¬ 
tween the two run states have to be attributed mostly to 
changes in bundle conformations leading to changes in 
the rotational drag . Large variations in could re¬ 
sult from having one or more flagella rotating out of the 
bundle. From <§ and we obtain that the propelling 
force should be directly proportional to the measured 
body rotation frequency Fthrust = B^ /Since 
does not change upon bundle reversal, any change in 
the ratio Fthrust /has to be attributed to different bun¬ 
dle conformations in the two states leading to a different 
B^I ratio. Fig. [^a) reports a plot of Fthrust versus 
for all cells that undergo at least one reversal. Solid lines 
connect points that refer to the same cell. If the bun¬ 
dle conformation remained the same upon reversal, for 
each cell, the two points should lay on a line with slope 
D^B^/D^ passing through the origin. This is clearly 
not the case for the vast majority of cells. From data 
shown in Figl^we get A[Bf/Di] = A [Fthrust= 0.47 
indicating again a substantial variation between bundle 
conformations in the two states. The solid black line in 
Fig. [^a) represents the best fit through the points whose 
slope (0.0026 ± 0.0001 pN s/rad) provides the average 
value: 



where we have assumed that body size and flagel¬ 
lar bundle conformation are statistically independent. 
The average value of (F^) can be estimated approx¬ 
imating the cell body with a prolate ellipsoid having 
the cell length I as the major axis and the cell thick¬ 
ness a as the average minor axis, (where 

fi is the water viscosity). The average cell length and 
its corresponding standard error are obtained by 3D re¬ 
constructions ((^) = 3.12 ± 0.06 jam) while the width 
is independently measured on fluorescently labeled cells 
((a) = 0.40 ± 0.02 jam). As a result we obtain an av¬ 
erage value for (D^) = 0.0063 ± 0.0006 pN /am s/rad. 
Using ^ we find {B^/D^) = 0.41 ± 0.04 jam Our 
value is significantly lower then that found in [T0| , where 
the authors obtained for the same ratio 1.13 /im“^. Re- 
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cent studies [21] have shown that resistive force theories 
usually result in an overestimation of the thrust of 
helical propellers while Lighthill’s slender body theory 
[22] provides a much more accurate result. In particular, 
using average bundle parameters from [2] (length 7.3 /im, 
pitch 2.3 yum, diameter 0.36 jam) and assuming a bundle 
radius of 0.02 jam we get from slender body theory a ratio 
/D^ =0.38 that is in very good agreement with our 
measured value. 
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FIG. 4. (a) Thrust idhrust versus body angular frequency 

for all the bacteria undergoing bundle reversals. Different 
bacteria are shown in different colours. Solid lines join the two 
run states for the same cell. Error bars represent ± the stan¬ 
dard deviation of idhrust and between different runs, (b) 
Histogram of the normalized force difference between the two 
run states, (c) Histogram of the normalized torque difference 
between the two run states. 


CONCLUSIONS 


We combined optical traps and 3D holographic mi¬ 
croscopy to investigate asymmetric features that arise in 
E. coll swimming as a consequence of flagellar bundling 
at either pole of the cell. We found that two well defined 
running states exist. These are characterised by large 
differences in propelling force, total applied torque and 
flagellar conformation. We also provide a direct mea¬ 
surement of the ratio B^ /which confirms predictions 
from slender body theory. Our technique can be easily 
extended to multiple cell systems allowing to study hy¬ 
drodynamic couplings between individual swimming cells 
in well defined and reproducible configurations. 

The research leading to these results has received 
funding from the European Research Council under 
the European Union’s Seventh Framework Programme 
(FP7/2007-2013) / ERC grant agreement n° 307940. We 
also acknowledge funding from MIUR-FIRB project No. 
RBFR08WDBE. 


APPENDIX: TRAP CALIBRATION 





FIG. 5. Probability distribution (a-c) of x, z, cos(0) of a 
trapped non-motile E. coli cell. The corresponding potentials, 
measured in units of ksT, are plotted as blue dots (d-f). The 
potentials are fitted with polynomials of respectively second, 
third, and first order that are plotted as blue lines. 


Trap calibration is performed individually on each cell 
after a short exposure to 532 nm laser light that stops 
motility. The Brownian fluctuations of position x^y^z 
and polar angle 0 are recorded for 5 minutes in an op¬ 
tical trap of reduced intensity such that fluctuations in 
^ have a comparable amplitude to those observed with 
motile cells. For small fluctuations the optical energy is 
separable in the sum of independent contributions: 


U{x,y,z,e) ~ UFx) + Uy{y) + U^z) + U0{0)) (7) 

Thermal fluctuation will be then distributed according to 
Boltzmann statistics leading to a factorized probability 
density: 

P(x, y, z, 9) oc (8) 

with p = IjkBT. This approximation is confirmed by the 
absence of correlation between the coordinates (x, y, z, 0). 
The individual probability densities and the correspond¬ 
ing potentials are reported in Fig. We are interested 
here in converting 2 : displacements into optical forces. To 
this aim we fit Uz{z) with a 3^^ order polynomial: 

U{Z) = ^kz^ + (g) 

We find that k'/k is always of the order of 10“^ /im“^. 
In our experiments the laser power is adjusted to the 
minimum value that is required for stable trapping of 
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swimming cells. In this situation the cell body moves 
away from the trap’s equilibrium position by 1.2 fim at 
most. Therefore, the third order term only affects force 
measurements at the level of a few percent that is below 
other sources of error. 
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